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ABSTRACT 



Aims. Recent studies have detected linear polarization in L dwarfs in the optical I band. Theoretical models have been developed to 
explain this polarization. These models predict higher polarization at shorter wavelengths. We discuss the polarization in the R and I 
band of 4 ultra cool dwarfs. 

Methods. We report linear polarization measurements of 4 ultra cool dwarfs in the R and I bands using the Intermediate dispersion 
Spectrograph and Imaging System (ISIS) mounted on the 4.2m William Herschel Telescope (WHT). 

Results. As predicted by theoretical models, we find a higher degree of polarization in the R band when compared to polarization in 
the I band for 3/4 of these ultra cool dwarfs. This suggests that dust scattering asymmetry is caused by oblateness . We also show how 
these measurements fit the theoretical models. A case for variability of linear polarization is found, which suggests the presence of 
randomly distributed dust clouds. We also discuss one case for the presence of a cold debris disk. 

Key words, polarization - atmospheric effects - brown dwarfs - low mass stars 



1. Introduction 

A large number of ultra cool dwarfs have been detected in 
the last decade, and our understanding of these faint ob- 
jects has kept improving. One of the challenging and fun- 
damental aspects in the study of these objects is to under- 
stand the properties and distribution of condensate dust in 
the atmosphere. Observations of L dwarfs with effective tem- 
peratures of 1400-2200 K have led to the investigation of 



• dust condensates in their atmospheres ( Kirkpatrick et. al. 1999 
[ |Tsuji et al. 19 96'). Because of complete gravitational settling, 
■ grains are expected to condense beyond the visible atmo- 
sphere for objects with effective temperatures below 1400 
K (T-Dwarfs - lAUard et al. 2"00n IChabrier et al. 2000l l. At 
higher effective temperatures (1400-2200 K), grains can be 
present in the visible atmosphere because of incomplete grav- 
itational settling (Burrows & Sharp 1999( IBurrows et al. 20011 



rotation of L dwarfs will induce the shape of their photo- 
sphere into the form of an oblate ellipsoid.! Basri et al. 20001) 
and this nonsphericity will lead to the incomplete cancel- 
lation of the polarization from different areas of the stel- 
lar surface ( Sengupta & Krishan 200 1[ ). This prediction was 
first confirmed by the detection of linear polarization at 768 
nm from a few L dwarfs by Menard et al. 2002 Recently, 
|Zapatero Osorio et al. 2005] have reported R and I band detec- 
tion of linear polarization from several L dwarfs. Since polar- 
ization in the optical is unlikely to be due to Zeeman splitting 
of atomic or molecular lines or by synchrotron radiation, the 
observed polarization can be explained by single dust scatter- 
ing in a rotationally induced oblate atmosphere ( jSengupta 2003[ 



Ackerman & Marley 2001[ lAllard et aL2bOH |Tsuji et al. 2004 



Cooper et al. 2003 Helling et al. 2008 ). Recent discoveries of 



blue L dwarfs and L-T transition type dwarfs (as identified 



in|Knapp et al. 2004 [IChiu et al. 2006l|Tsuji & Nakajima 20031 



have brought forth models which could explain this phenomenon 



(eg. Burrows, Sudarsky & Hubeny 2006 Knapp et al. 2004 i by 
mechanisms which involve dust settling. It would be very impor- 
tant to validate these mechanisms. 

Linear polarization could be a very useful tool in under- 
standing the observationally poorly constrained dust properties 
in the atmospheres of L dwarfs. The possibility of detecting po- 
larization at optical wavelengths from grains in the atmospheres 
of L dwarfs was first raised by [Sengupta & Krishan 200 1| Fast 



Sengupta 8c Kwok 2005| ) or it could be due to la rge and ran- 
domly distnbuted dust clouds JMenard et al. 2002] i. 



In this paper, we report polarization measurements of 3 L 
dwarfs (L0-L5) and one M9.5 dwarf with WHT/ISIS in both I 
and R bands . We also discuss our results comparing them with 
the recently published results of lGoldman et al. 20091 



Our measurements show the general trend that polarization 
is higher in the R band than the one in the I band. This trend 
strongly supports the presence of dust in the atmosphere of L 
dwarfs as it is very unlikely that any other mechanisms (such as 
the presence of magnetic field) can explain this observation at 
optical wavelengths (Mena rd et al. 2002] l. We also discuss how 
the theoretical models (see Sect. 3) successfully fit our measured 
data. 
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2. Target selection and observations 

Four very nearby (7pc -15pc) ultra cool dwarfs were selected 
(SpT M9.5 - L5) as they are among the brightest and nearest 
ultra cool dwarfs with no known infrared excess and no evidence 
of multiplicity(see Table 1). These selection criteria ensure that 
the targets are bright enough sources in the R band to get high 
S/N and to avoid other than intrinsic sources of polarization 
such as circumstellar disks or multiplicity. For calibration, one 
polarized (Cyg OB2 A Whittet et.al. 1992) and one unpolarized 
standard star were observed at two different times during the 
night. All the objects were observed in a way that they insured 
they were acquired at the same position on the detector (5 
pixel box). This procedure minimized contamination caused by 
instrumental polarization within the detector and variations in 
the optical path. 

The polarimetric observations were obtained using the 
Intermediate dispersion Spectrograph and Imaging System 
(ISIS) which is mounted at the Cassegrain focus of the 4.2m 
William Herschel Telescope (located in La Palma, Canary 
Islands, Spain). ISIS in polarization mode is a modulation 
polarimeter with a double-beam analyzer (the calcite plate) and 
a rotating halfwave plate modulator ISIS is equipped with two 
detectors; a blue-sensitive EEV12 (4096x2048 pixels) and a 
red-sensitive RED-i- (4096x2048) detector. In our program, we 
have used the RED+ detector. 

Images were obtained using Bessel R- and I- filters centered 
on 641 and 812 nm, respectively, on June 18, 2006 (UT Date). 
The night was photometric with stable average seeing of 1.0 arc- 
sec . 

The raw images were bias-subtracted and flat-fielded before 
performing aperture photometry. The flat-field images were ob- 
tained with the polarimeter optics. Fluxes were obtained for 0.8, 
1.0, 1.2, 1.5, 2.0 times the average FWHM for each object. The 
best aperture was chosen to be 1 .5 times FWHM based on mini- 
mum photon contribution of nearby sources, variable sky contri- 
bution, and maximum signal-to-noise ratio of the measurements. 
The average FWHM of all images was 4.0 pixels which corre- 
sponds to 1.0 arc sec. We have only one set of measurements for 
each object. Therefore, we have estimated the uncertainty in the 
degree of polarization from various apertures (a similar method 
was used by Zapatero Osorio et al. 2005 for some of the ob- 
jects). There was no significant instrumental polarization found 
as the unpolarized standard measured D(p)=0.086% + 0.002 . 

3. Theoretical modeling of polarization 

Polarization is a measure of anisotropy in the radiation field and 
is caused by either scattering or is due to the presence of mag- 
netic field. The state of the polarization of light is described by 
the Stokes parameters, /, Q, U and V. The parameter / is the total 
scalar specific intensity of radiation. It is the complete flux of 
radiant energy inside the unit intervals of frequency, time, solid 
angle, and area perpendicular to the flux. This flux includes all 
radiation independently of polarization. Polarization is described 
by the parameters Q, U, V. These parameters are proportional to 
the scalar specific intensity and have the same dimension. Q and 
U represent the linearly polarized component, and V represents 
the circularly polarized component. For linear polarization, V-0 
and the degree of polarization is given as /? = VG^ + C/^//. If 
we consider axial symmetry, then U - Q and in that case we de- 
fine the degree of polarization p - -QII. The sign convention 



is such lliat if /? > 0, the light is polarized perpendicular to the 
scattering plane, and if p < 0, the light is polarized parallel to 
the scattering plane. For an unresolved stellar object, the Stokes 
parameters are integrated over the stellar disk. 

From our obtained images, the degree of polarization and the 
polarization angles are calculated using the following equations: 
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o(45°)M45°) 
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where 



o,e ordinary and extraordinary fluxes respectively of 
the dual images of the program source on a single 
frame 

p degree of polarization 

6 polarization angle 

/ the total scalar specific intensity of radiation 

Q/I,U/I normaUzed Stokes parameters 



As pointed out by IMenard et al. 20021 the observed lin- 
ear polarization in the optical cannot be due to magnetic 
field, and scattering remains the most viable mechanism for 
yielding the detected linear polarization. Polarimetric obser- 
vation at the R and I bands by 'Zapatero Osorio et al. 2005 
shows that polarization decreases significantly with the in- 
crease in wavelength, which strongly supports this argument 
( Sengupta & Kwok 2005, ). In the present investigation, we re- 



port detection of polarization at both R and I bands which shows 
the same wavelength dependency and hence strengthens the case 
for scattering polarization. If the dust density is low, then the sin- 
gle scattering approximation is reasonable for the region where 
the dust optical depth < 1 because scattering by atoms 
and molecules does not significantly contribute to polarization. 
However, in the optical region is much higher than 1, and 
therefore multiple scattering is the appropriate process. Multiple 
scattering reduces the amount of polarization significantly be- 
cause the planes of scattering are oriented randomly and can- 
cel each others' contributions. In the present paper, however , 
we adopt a simple single scattering polarization model because 
the main aim is to explain the wavelength dependency of the 
observed polarization. Nevertheless, we find that a single dust 
scattering model in rotation-induced oblate L dwarfs can fit the 
observed data well within the permissible range of parameter 
space of dust composition, size, geometrical distribution and ro- 
tation period. A detailed modeling with multiple scattering is in 
progress. 

The simple theoretical model adopted 

§ lain the observed po larization is described 
engupta & Kwok 2005 
larization integrated over the stellar disk is given by : 



here to ex- 
in details in 
At an edge-on view, the degree of po- 



P(A): 



(A2- 1)^2] 1/2 



dfi\(5) 



In the above expression, P is the total pressure (gas plus dust), 
P,op and Phot are the pressures at the top (deck) and the bottom 
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(base) of the dust layer, n(P) is the grain number density at dif- 
ferent pressure heights, p{P) is the total density at any pressure 
level, and A is the ratio of the equatorial radius to the polar ra- 
dius. Since the dust density is small compared to the gas density, 
p{P) is the effective gas density and is equal to nRT/P where 
n is the mean molecular weight and R is the gas constant. In 
the above equation, Pi(i-i) is the Lengendre polynomial where 
ju - cos 6, 6 being the scattering angle and ¥{1 is given by 



¥lm 

where 

a(l, m) 



Pf(cos0)£/(cos6l). 



(21 + - m) 



4-n{l + m) 



1/2 



(6) 



(7) 



and P'" is the associated Legendre function of the first kind, ii 
and 12 are the scattering functions given bv lvan de Hulst 19571 

The vertical dust distribution and the location of the cloud 
base and deck in the atmosphere are calculated based on the one 



dimensional heterogeneous cloud model of Cooper et al. 2003 



This model assumes chemical equilibrium throughout the atmo- 
sphere and uniform density distribution across the surface of an 
object at each given pressure and temperature. The number den- 
sity of cloud particles in this model is given by 



(8) 



where p is the mass density of the surrounding gas, a is the cloud 
particle radius, pd is the mass density of the dust condensates, fi 
and fid are the mean molecular weight of atmospheric gas and 
condensates respectively. The condensate mixing number ratio 
iqc) is given as 



~ q below ' 



(9) 



for heterogeneously condensing clouds. In the above equation, 
Ibeiow is the fraction of condensible vapor just below the cloud 
base, PcJ is the pressure at the condensation point. 

A log-normal size distribution is adopted for the spheri- 
cal grains ( Ackerman & Marley 2001 1. The formation of dust 
makes it a prohibitive task to develop a fully consistent atmo- 
spheric model for ultra-cool dwarfs. This is mainly because of 
the fact that the presence of dust clouds affects the radiative 
equilibrium of the upper atmosphere and hence alters the T- 
P profile from that of a cloud-free atmosphere. On the other 
hand, the T-P profile dictates the position and the chemical equi- 
librium of condensates. Allard et al. (2001) presented atmo- 
spheric models for two of the limiting cases, e. g., one with 
inefficient gravitational settling wherein the dust is distributed 
according to chemical equilibrium predictions (AMES-dusty) 
and another with efficient gravitational settling in which sit- 
uation dust has no effect on the thermal structure (AMES- 
cond). [Tsuji et al. 2004| have proposed a Unified Cloudy Model 
(UCM) in which the segregation of dust from ffie gaseous mix- 
ture takes place in all the ultra-cool dwarfs and at about the 
same critical temperature. [Ackerman & Marl ey 200 1| treat the 
upward convective mixing of a gas, its condensation and the sed- 
imentation of the condensate through the atmosphere of the ob- 
ject while |Woitke & Helhng 2004| consider an ensamble of dust 



grains falUng downwards from the top of the atmosphere. A de- 
tailed comparison of different atmospheric models of L dwarfs 



is presented in Helling et al. 2008 



The oblateness of a rotating object has been discussed by 
IChandrasekhar 2003] in the context of poly tropic gas configura- 
tion under hydrostatic equilibrium. For a slow rotator, the rela- 
tionship for the oblateness / of a stable polytropic gas configu- 
ration under hydrostatic equihbrium is given by 



2 Q?Rj 
f = -c -, 

■' 3 GM ' 



(10) 



where M is the total mass, R^ is the equatorial radius, and Q 
is the angular velocity of the object which is related to the lin- 
ear velocity V = Q. x R,,. C is a constant whose value depends 
on the polytropic index. For a polytropic index of n = 1.0, 
C - 1.1399, which is appropriate for Jupiter (IHubbard 19841 1. 
For non-relativistic completely degenerate gas, n - 1.5 and 
C = 0.9669. 

The effective temperature of the L dwarfs of different spec- 
tral type is determined by adopting a sixth order polynomial 
fit given bv Golim owski et al . 2004' which is based on bolomet- 
ric luminosities. The T^ff calibration of Golimowski et al. 2004] 
agrees well in the interval L3-L8, but there are significant differ- 
ences in earlier types. In our calculations for the degree of po- 
larization, the effective temperature T^ff is used and hence the 
degree of polarization should be considered strictly as a function 
of Teff rather than of spectral type. The mass and radius of the L 
dwarfs of different spectral types are estimated by adopting the 
empirical relationship given by |Marley et al. 1996 



4. Results 

We find a trend (more data are required to confirm our theory) of 
higher polarization in the R band when compared to the I band. 
This wavelength dependency strongly supports the argument by 
|Sengupta & Kwok 2005| that the polarization arises due to scat- 
tering and not because of magnetic field. In dust scattering as 
described by Mie theory, the amount of polarization depends on 
the ratio of the grain radius to the wavelength. For the same kind 
of dust species, the polarization usually peaks when the ratio 
is one. As a consequence, the increase in polarization with the 
decrease in wavelength implies the presence of sub-micron size 
grains in the photosphere of ffie L dwarfs. We present our mea- 
surements in Table 2, and our model fit in Figure 1 . One of the L 
dwarfs from Zapatero Osorio et al."2005l (2MASSW J1507476- 
162738) shows a null polarization in our measurements in the I 
band, whereas Zapatero Osorio et al. 2005 present a higher po- 
larization (1.36% ±0.30) for ffie same object. This suggests vari- 
ability in linear polarization which in turn suggests atmospheric 
activities like dynamical variations of the cloud cover. 

We also find relatively high polarization in the LO dwarf. 
Additionally, from iSchmidt et al. 20071 ffie Ha equivalent width 
is the highest for the LO dwarf among the objects from our 
sample. This could be an indirect evidence of a disk around 
this ultra cool dwarf. We therefore searched the Spitzer public 
archive for mid-IR data. 2MASS J17312974-1-2721233 has been 
observed with IRAC and IRS in the course of program 3136 
(P.I. Cruz), and we retrieved the pipeline processed data. We 
extracted the IRAC photometry using standard PSF photometry 
procedures within the Interactive Data Language. Uncertainties 
were estimated from the Poisson noise weighted by the cov- 
erage maps of the mosaics. Table |4] gives a summary of the 
photometry. Figure |2] shows the spectral energy distribution 



4 

Table 1. Target list 
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Object 


RA 


DEC 


SpT 


R 


I 


J 


d 


Ha 


ref 










(mag) 


(mag) 


(mag) 


(pc) 


EW(A) 




2MASSW J1438082+640836 


14 38 08.2 


+64 08 36 


M9.5 




16.8 


12.92 


7.2 +0.5 


-2.36 


1 


2MASSW J1507476-162738 


15 07 47.69 


-16 27 38.6 


L5 


18.9 


16.5 


12.8 


7.33 ± 0.03 


-5.01 


3,4,5 


2MASS J17312974+2721233 


17 31 29.74 


+27 21 23.3 


LO 






12.09 


11.80 + 0.70 


-5.98 


1,2 


2MASSI J1807159+501531 


18 07 15.93 


+50 15 31.6 


L1.5 






12.96 


14.6 ± 1.0 


-1.38 


1,2,3 



Note 1. ref- l.lSchmidt et al. 20071 2. IJameson et al. 20071 3. ICruz et al. 2003] 4. IReid et al. 20001 5. |Knapp et al. 2004 | 
Ho- EW from lSchmidt et al. 20071 



(SED) of the source and a LO comparison object from the 
literature (2MASS J1204+3212, Patten et al., 2006). 2MASS 
J 173 12974+272 1233 does not show any significant mid-IR ex- 
cess up to 15 //m. The presence of a young circumstellar disc can 
therefore be ruled out at a high level of confidence. In the current 
state of the data, we cannot rule out the presence of a cold debris 
disc, as it would produce an excess at longer wavelengths. 

The surface gravity of L dwarfs older than approximately a 
few hundred million years varies from g ^ 10^ to 3 X 10^ cms . 
Evolutionary models by Chabrier et al. 20001 show younger L 
dwarfs to have a surface gravity smaller than 10^ cms^^. Our 
best fit model parameters are presented in Table 3 (best fit 
obtained by visual inspection). For the L dwarfs 2MASSW 
J1438082+640836, we find the best fit with log(g)=5.0, and ro- 
tational velocity V=18.7 km s"' when n - 1.0, and V=20.5 km 
s"' when n - 1.5. For a fixed rotational velocity, the oblate- 
ness is lower when the polytropic index n is higher. As a result, 
higher rotational velocity is needed to fit the observed data if n = 
1.5. The projected rotational velocity of the L dwarf 2MASSW 
J1507476-162738 is measured by Bailer-Jones 2004 while the 
projected rotati onal velocity of 2MAS SW J1807159+501531 
is estimated by IReiners & Basri 20081 In the absence of any 
knowledge on the projection angle, we consider the minimum 
rotational velocity of these objects as V=27.2 km s"' and 
V=76 km s We find the best fit for the observed polar- 
ization of 2MASSWW J1507476-162738 with Log(g)=5.255 
for n = 1.0 and Log(g)=5.230 for n = 1.5 For 2MASSW 
J1807159+501531 the best fit is obtained with log(g)=5.41, 
when n - 1.0, and log(g)=5.37 when n - 1.5. Note that the 
oblateness decreases with the increase in surface gravity. For all 
the cases, the observed polarization profiles can be fitted with 
sub-micron size grains, and the mean size of grains that are re- 
quired to fit the observation is consistent with the recent theo- 



pared to the I band as predicted by the theoretical models of 



Sengupta & Kwok 2005 



retical calculations of dust properties ( Woitke & Helling 2004 



Woitke & Helling 2003 1. Polarization measurements for one 
of the above three objects (2MASS W J1507476 -162738) 
were also recently published Goldma n et al. 20091 Our re- 
sults are consistent with the Goldma n et al. 2009| measurements 
within lo" error bars for this object. This work differs from 
IGoldman et al. 20091 in the aspect that our sample could re- 
produce the predictions of' Sengupta & Kwok 2005| whereas the 
jGoldman et al. 2009 sample does not. Both the studies (this pa- 
per and IGoldman et al. 2009) have sparse statistics, which dis- 
courages us to make any speculations on this apparent discrep- 
ancy. 

5. Conclusions 

1 . We report linear polarizaion measurements of 4 very nearby 
ultra cool dwarfs in the R and I bands. 

2. We find that there is a trend (3 out of 4) of a higher degree 
of polarization at shorter wavelengths (R band) when com- 



3. The LO dwarf 2MASS J17312974+2721233 is interest- 
ing because of its relatively high polarization and requires 
follow-up studies. 

4. We also fit theoretical models to predict the dust grain size 
and rotational velocities of three of the ultra cool dwarfs. 

5. We find evidence for variability in the linear polarization for 
(2MASSW J1507476-162738). This suggests atmospheric 
activities like dynamical variations of the cloud cover in this 
object. 
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Table 2. Polarization measurements 



5 



WDJcCl 


o pi 


riiier 


Exposure Tinie 




u 








(seconds) 


(%) 


(degrees) 






T 
1 


on 




UU.Z.H- 






R 


300 


0.54 ±0.097 


115 


2MASSW J1507476-162738 


L5 


I 


150 


0.0 ±0.036 








R 


300 


0.216 ±0.022 


100.26 


2MASS J17312974+2721233 


LO 


I 


60 


5.195 ± 0.854 


52.57 






R 


240 


0.666 ±0.169 


111.08 


2MASSIJ1807159+501531 


L1.5 


I 


60 


0.711 ±0.142 


46.67 






R 


120 


1.669 + 0.605 


101.2 



Note 2. The uncertainty in due to calibration errors is about 1.2 degrees 
Table 3. Model fit 



Object 


SpT 


n 


do 


log(g) 


V 








(micron) 


(cgs) 


(kms-i) 


2MASSW J1438082+640836 


M9.5 


1.0 


0.7 


5.0 


18.7 






1.5 


0.7 


5.0 


20.5 


2MASSW J1507476-162738 


L5 


1.0 


0.44 


5.255 


27.2 






1.5 


0.44 


5.230 


27.2 


2MASSI J1807159+501531 


L1.5 


1.0 


0.5 


5.41 


76.0 






1.5 


0.44 


5.37 


76.0 



Note 3. Bailer-Jones 2004 has measured the projected rotational velocity for 2MASSW J1507476-162738 as 27.2 km s"', while Reiner & Basri 
2008 have reported the projected velocity for 2MASSW J1807159+501531 to be 76 km s"'. 

Table 4. IRAC photometry of 2MASS J173 12974+2721233 



Wavelength 


Flux 


(micron) 


(mJy) 


3.6 


21.82±0.02 


4.5 


14.16±0.02 


5.8 


11.30+0.06 


8.0 


6.91±0.03 




Fig. 1. Best model fit of the observed data. The solid Unes represent the model with the polytropic index n=1.0 and dashed lines 
represent that with n=1.5. For other parameters see Table 3. 
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Fig. 2. Spectral energy distribution of 2MASS J173 12974+2721233 (dots). V-band photometry from the LSPM-North proper- 
motion catalog of nearby stars (Lepine et a., 2005). J,H and Ks photometry from 2MASS. IRAC and IRS photometry from this 
work. A comparison of LO dwarf (2MASS J1204+3212, Patten et al., 2006) scaled to the same J-band flux is overplotted (squares). 
The target does not show any significant excess compared to the field LO dwarf. 



